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Abstract. The emission of gravitational waves from a system of massive objects interact- 
ing on elliptical, hyperbolic and parabolic orbits is studied in the quadrupole approximation. 
Analytical expressions are then derived for the gravitational wave luminosity, the total en- 
ergy output and gravitational radiation amplitude. A crude estimate of the expected number 
of events towards peculiar targets (i.e. globular clusters) is also given. In particular, the rate 
of events per year is obtained for the dense stellar cluster at the Galactic Center. 
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1. Introduction 

Gravitational waves are generated by dynam- 
ical astrophysical events, and they are ex- 
pected to be strong enough to be detected 
when compact stars such as neutron stars (NS) 
or black holes (BH) are involved in such 
events. In particular, coalescing compact bi- 
naries are considered to be the most promis- 
ing sources of gravitational radiation that can 
be detected by the ground-based laser interfer- 
ometers. Advanced optical configurations capa- 
ble of reaching sensitivities slightly above and 
even below the so-called standard-quantum- 
limit for a free test-particle, have been de- 
signed for second and third generation GW de- 
tectors. A laser-interferometer space antenna 
(LISA) (1(T 4 ~ lQ- 2 Hz) might fly within the 
next decade. It is important in order to predict 
the accurate waveforms of GWs emitted by ex- 
treme mass-ratio binaries, which are among the 
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most promising sources for LISA. To this aim, 
searching for criteria to classify the ways in 
which sources collide is of fundamental impor- 
tance. A first rough criterion can be the classi- 
fication of stellar encounters in collisional as in 
the globular clusters and 'mcollisio nless as in 
the galaxies iJ.Binney et al. ] d!987h . A funda- 
mental parameter is the richness and the den- 
sity of the stellar system and so, obviously, we 
expect a large production of GWs in rich and 
dense systems. Systems with these features are 
the globular clusters and the galaxy centers. In 
particular, one can take into account the stars 
(early-type and late-type) which are around the 
Galactic Center, Sagittarius A* (SgrA*) which 
could be very interesting targets for the above 
mentioned ground-based and space-based de- 
tectors. In recent years, detailed information 
has been achieved for kinematics and dynam- 
ics of stars moving in the gravitational field 
of such a central object. The statistical prop- 
erties of spatial and kinematical distributions 
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are of particular interest. Considering a field of 
resolved stars whose proper motions are accu- 
rately known, one can classify orbital motions 
and deduce, in principle, the rate of produc- 
tion of GWs according to the different types of 
orbits. This work is motivated by a classifica- 
tion of orbits in accordance with the conditions 
of motion, we want to calculate the GW lumi- 
nosity for the different types of stellar encoun- 
te rs. A similar approach has been developed 
in IS. Capozziello and M. De Laurentis I (120081) 
but, in that case, only hyperbolic trajectories 
have been considered. In this report we in- 
vestigate the GW emission by binary systems 
considering bounded (circular or elliptical) and 
unbounded (parabolic or hyperbolic) orbits. 
We expect that gravitational waves are emit- 
ted with a "peculiar" signature related to the 
encounter-type: such a signature has to be a 
"burst" wave-form with a maximum in cor- 
respondence of the periastron distance. The 
problem of bremsstrahlung-\\ke gravitational 
wave emission has been studied in detail by 
Kovacs and Thorne by considering stars inter- 
acting on unbounded and bounded orbits. In 
this report, we face this problem discussing 
in detail the dynamics of such a phenomenon 
which could greatly improve the statistics of 
possible GW sources. 

2. Orbits in stellar encounters 

Let us take into account the Newtonian the- 
ory of orbits since stellar systems, also if at 
high densities and constituted by compact ob- 
jects, can be usually assumed in Newtonian 
regime. We give here a self-contained sum- 
mary of the well-known orbital types in order 
to achieve below a clear classification of the 
possib le GW emissions. We refer to the text 
books iJ.BinnevetaTl dl987l) : |L. Landau eTall 
(119731) for a detailed discussion. A mass m\ is 
moving in the gravitational potential <I> gener- 
ated by a second mass ni2- The vector radius 
and the polar angle depend on time as a con- 
sequence of the star motion, i.e. r = r(f) and 
<p = <p(f). the total energy and the angular mo- 
mentum, read out 
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respectively, where u = m ' m2 is the reduced 

r ' m\+m2 

mass of the system and y - Gm\mi. To solve 
the above differential equations we write 

dr dr d(p L dr L d 1 1 \ 
dt d<p dt fir 2 dip fi d(f>\r J 

and defining, as standard, the auxiliary variable 
u = 1/r, Eq. (HJ takes the form 
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where u' = du/d<p and we have divided by 
L 2 /2fi. Differentiating with respect to <p, we get 



u' (u" + u - ^) = 



(5) 



hence either u' = 0, corresponding to the cir- 
cular motion, or 



u + u 



W 
L 2 



which has the solution 
yu 

u — — + C cos ((f> + a) 
L 1 

or, reverting the variable, 



L 2 



C cos (<p + a) 



(6) 
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(8) 



which is the canonical form of conic sections 
in polar coordinates. The constant C and a are 
two integration constants of the second order 
differential equation ©. The solution ([8]l must 
satisfy the first order differential equation 
Substituting © into (0)1 we find, after a little 
algebra, 



r 2 _ 2fiE lyiiY 
~ L 2 \L 2 ) 



(9) 



and we get C 2 > 0. This implies the 
four kinds of orbits given in Tabl e I (se e 
dS. Capozziello and M. De Laurentis! l2008l) ). 
Circular motion correspond to 
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C = E = E m j„ circular orbits 

< |C| < ff~ £„„■„ < E < elliptic orbits 

\C\ = jj E = parabolic orbits 

\C\ > j| E > 0, hyperbolic orbits 



well known that the amplitude of GWs can be 
evaluated by 

hj k (t,R)=—& k , (13) 

R being the distance between the source and 
the observer and { j, k] = 1,2, where is the 
quadrupole mass tensor 

a7 = 2 mfl(3x « x °"^ 2) ' (14) 

a 

. Here G being the Newton constant, r a the 
modulus of the vector radius of the a - th par- 
ticle and the sum running over all masses m a 
in the system. We now derive the GW ampli- 
tude in relation to the orbital shape of the bi- 
nary systems. As an example, the amplitude of 
gravitational wave is sketched in Fig. Q] for a 
stellar encounter close to the Galactic Center. 
The adopted initial parameters are typical of a 
close impact and are assumed to be b — 1 AU 
and v'o = 200 Km.? -1 , respectively. Here, we 
have fixed M\ = M 2 = 1.4M . The impact pa- 
rameter is defined as L = bv where L is the 
angular momentum and v the incoming veloc- 
ity. We have chosen a typical velocity of a star 
in the galaxy and we are considering, essen- 
tially, compact objects with masses compara- 
ble to the Chandrasekhar limit (~ 1 .4M ). This 
choice is motivated by the fact that ground- 
based experiments like VIRGO or LIGO ex- 
pect to detect typical GW emissions from the 
dynamics of these objects or from binary sys- 
tems composed by them. 



Table 1. Orbits in Newtonian regime classified 
by the approaching energy. 
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Elliptic motion 
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where e — J is the eccentricity of the el- 
lipse and / is semi-latus-rectum of the ellipse. 
Parabolic and Hyperpolic motion correspond 
to 

l + ecos^>0 (12) 

This means cos <p > -1, i.e. e (-n, n) and the 
trajectory is not closed any more. For <p — > ±n, 
we have r — > oo. The curve, with e = 1, is 
a parabola. For e > 1, the allowed interval of 
polar angles is smaller than (f> € {-n, n), and 
the trajectory is a hyperbola. Such trajectories 
correspond to non-returning objects. 



3. Gravitational wave 

At this point, considering the orbit equa- 
tions, we want to classify the gravitational 
radiation for the different stellar encounters 
dC. W Misner et al. I fl973h . Direct signatures 
of gravitational radiation are its amplitude and 
its wave-form. In other words, the identifica- 
tion of a GW signal is strictly related to the ac- 
curate selection of the shape of wave-forms by 
interferometers or any possible detection tool. 
Such an achievement could give information 
on the nature of the GW source, on the propa- 
gating medium, and , in principle, on the gravi- 
tational theory producing such a radiation. It is 



4. Rate and event number 
estimations 

An important remark is due at this point. 
A galaxy is a self-gravitating collisionless 
system where stellar impacts are very rare 
dJ.Binnev et al. |[l987l) . From the GW emission 
point of view, close orbital encounters, colli- 
sions and tidal interactions should be dealt on 
average if we want to investigate the gravita- 
tional radiation in a dense stellar system as we 
are going to do now. 

Let us give now an estimate of the stellar 
encounter rate producing GWs in some inter- 
esting astrophysical conditions like a typical 
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Fig. 1. The gravitational wave-forms from el- 
liptical orbits shown as function of the polar 
angle 4>. We have fixed M\ = M 2 = 1 .4M Q . M 2 
is considered at rest while Mi is moving with 
initial velocity vo = 200 Kms -1 and an impact 
parameter b — 1 AU. The distance of the GW 
source is assumed to be R — 8 kpc and the ec- 
centricity is e = 0.2, 0.5, 0.7. 



Fig. 2. The gravitational wave-forms for a 
parabolic encounter as a function of the polar 
angle ip. As above, M\ = M 2 = \ AM @ and M 2 
is considered at rest. Mi is moving with ini- 
tial velocity Vo = 200 Km^ 1 with an impact 
parameter b — 1 AU. The distance of the GW 
source is assumed at R = 8 kpc. The eccentric- 
ity is e = 1 . 



globular cluster or towards the Galactic Center 
after we have discussed above the features dis- 
tinguishing the various types of stellar encoun- 
ters. Up to now, we have approximated stars 
as point masses. However, in dense regions 
of stellar systems, a star can pass so close to 
another that they raise tidal forces which dis- 
sipate their relative orbital kinetic energy. In 
some cases, the loss of energy can be so large 
that stars form binary or multiple systems; in 
other cases, the stars collide and coalesce into 
a single star; finally stars can exchange grav- 
itational interaction in non-returning encoun- 
ters. To investigate and parameterize all these 
effects, we have to compute the collision time 
t co ii, where 1 / t co u is the collision rate, that is, 
the average number of physical collisions that 
a given star suffers per unit time. For the sake 
of simplicity, we restrict to stellar clusters in 
which all stars have the same mass m. Let us 
consider an encounter with initial relative ve- 
locity vo and impact parameter b. The angular 
momentum per unit mass of the reduced parti- 
cle is L = bvo. At the distance of closest ap- 
proach, which we denote by r co #, the radial ve- 
locity must be zero, and hence the angular mo- 



mentum is L = r co uVmax, where v ma x is the rel- 
ative speed at r co u. From the energy equation 
0, we have 



b 2 = d« + 



4Gmr coU 



(15) 



If we set r c „u equal to the sum of the radii of 
two stars, then a collision will occur if and only 
if the impact parameter is less than the value of 
b, as determined by Eq.(fT3Tl. 

Let /(v fl )af 3 v fl be the number of stars per 
unit volume with velocities in the range \ a + 
d 3 \ a . The number of encounters per unit time 
with impact parameter less than b which are 
suffered by a given star is just /(v fl )c/ 3 v fl times 
the volume of the annulus with radius b and 
length vq, that is, 



f(\ a )nb 2 v d :i \ a 



(16) 



where vo = |v— v fl | and v is the velocity 
of the considered star. The quantity in Eq.dToli 
is equal to 1 /t co u for a star with velocity v: to 
obtain the mean value of l/t co u, we average 
over v by multiplying ( fTol i by /(v)/v, where 
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Fig. 3. The gravitational wave-forms for hyper- 
bolic encounters as function of the polar an- 
gle (p. As above, we have fixed M\ - M2 
I AM . M2 is considered at rest while Mi is 
moving with initial velocity Vo = 200 Kmf 1 
and an impact parameter b — 1 AU. The dis- 
tance of the source is assumed at R = 8 kpc. 
The eccentricity is assumed with the values 
e= 1.2,1.5,1.7. 

v = J f(\)(P\ is the number density of stars 
and the integration is over cPy. Thus 

_L = V f e -(v 2 +v2)/2^ 
t co ii 87T 2 CT 6 J 

/ 4Gmr co u\ 3 3 
r OT B |v - v fl | + d 3 \d\ a (17) 



We now replace the variable v fl by V = 
v - v„. The argument of the exponential is then 

la 2 , and if we replace the 



variable v by v„„ = v — — V (the center of mass 
velocity), then we have 



1 _ v 
t co ii ~ 8n 2 a- 6 

,-(v 2 c ,„+V 2 )/2o- 2 



4-Gmr coU 
r C oiiV + |rfV.(18) 



The integral over \ cm is given by 



3/ V . 



(19) 



Thus 
1 

tcoii 



-V^tr 2 

2cr 3 

(r 2 „V 3 +4GmVr C0 „)</V 



(20) 



The integrals can be easily calculated and 
then we find 



1 a r 2 
= 4 yTtvcnr 

t co u 



4 y/nvGmr c 



■oil 



(21) 



The first term of this result can be derived 
from the kinetic theory. The rate of interac- 
tion is vE (V), where S is the cross-section and 
(V) is the mean relative speed. Substituting 
S = 7ir 2 „;i and (V) = Act I yfn (which is appro- 
priate for a Maxwellian distribution whit dis- 
persion cr) we recover the first term of (f2TT> . 
The second term represents the enhancement in 
the collision rate by gravitational focusing, that 
is, the deflection of trajectories by the gravita- 
tional attraction of the two stars. 

If r* is the stellar radius, we may set r a ,u = 
2r». It is convenient to introduce the escape 



speed from stellar surface, 
to rewrite Eq. (l2TT) as 



1 



2Gm 



and 



r = — = 16 y/nvcrr 2 1 + —V| = 
tcoii \ 4o- 2 / 

16V7rvo7- 2 (l +0), 



where 







vj _ Gm 
4a 2 ~ 2o- 2 r« 



(22) 



(23) 



is the Safronov number dJ.Binnev et al. I 

Il987h . In evaluating the rate, we are consider- 
ing only those encounters producing gravita- 
tional waves, for example, in the LISA range, 
i.e. between 10~ 4 and 10~ 2 Hz . Numerically, 
we have 



r * 5.5 x 

3 



10-10 f v )(-ZL) 

UOkms-'At/A 2 / 



—J yrS 



© « 1 



(24) 
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r. 5 .5xicW^-) 2 f v \ 

\10 5 M o / UOkms- 1 / 

x(^)(^) J --' 9 »> < 25 > 

If » 1, the energy dissipated exceeds 
the relative kinetic energy of the colliding stars, 
and the stars coalesce into a single star. This 
new star may, in turn, collide and merge with 
other stars, thereby becoming very massive. As 
its mass increases, the collision time is shorten 
and then there may be runaway coalescence 
leading to the formation of a few supermas- 
sive objects per clusters. If << 1, much of 
the mass in the colliding stars may be liberated 
and forming new stars or a single supermassive 
objects. 

Note that when we have the effects of 
quasi-collisions in an encounter of two stars in 
which the minimum separation is several stel- 
lar radii, violent tides will raise on the surface 
of each star. The energy that excites the tides 
comes from the relative kinetic energy of the 
stars. This effect is important for © >> 1 since 
the loss of small amount of kinetic energy may 
leave the two stars with negative total energy, 
that is, as a bounded binary system. Successive 
peri-center passages will dissipates more en- 
ergy by GW radiation, until the binary orbit is 
nearly circular with a negligible or null GW ra- 
diation emission. 

Let us apply these considerations to the 
Galactic Center which can be modelled as 
a system of several compact stellar clusters, 
some of them similar to very compact globu- 
lar clusters with high emission in X-rays. 

For a typical compact stellar cluster 
around the Galactic Center, the expected event 
rate is of the order of 2 x 10 9 yrs 1 which 
may be increased at least by a factor 100 if 
one considers the number of globular clusters 
in the whole Galaxy eventually passing nearby 
the Galactic Center. If the compact stellar clus- 
ter at the Galactic Center is taken into account 
and assuming the total mass M 3 x 10 6 
M , the velocity dispersion cr 150 km s -1 
and the radius of the object R 10 pc (where 
= 4.3), one expects to have 10~ 5 open or- 
bit encounters per year. On the other hand, if 
a cluster with total mass M 10 6 M , cr 



150km s _1 and/? 0.1 pc is considered, an 
event rate number of the order of unity per year 
is obtained. These values could be realistically 
achieved by data coming from the forthcoming 
space interferometer LISA. As a secondary ef- 
fect, the above wave-forms could constitute the 
"signature" to classify the different stellar en- 
counters thanks to the differences of the shapes 
(see the above figures). 

5. Concluding Remarks 

We have analyzed the gravitational wave 
emission coming from stellar encounters in 
Newtonian regime and in quadrupole approx- 
imation. In particular, we have taken into ac- 
count the expected luminosity and the strain 
amplitude of gravitational radiation produced 
in tight impacts where two massive objects of 
1 AM Q closely interact at an impact distance of 
IAU. Due to the high probability of such en- 
counters inside rich stellar fields (e.g. globular 
clusters, bulges of galaxies and so on), the pre- 
sented approach could highly contribute to en- 
large the classes of gravitational wave sources 
(in particular, of dynamical phenomena capa- 
ble of producing gravitational waves). In par- 
ticular, a detailed theory of stellar orbits could 
improve the statistic of possible sources. 
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